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ABSTRACT 

We have determined accurate values of the product of the mass-loss rate and the ion fraction 
of P 4+ , Mq( P 4+ ), for a sample of 40 Galactic O-type stars by fitting stellar-wind profiles to 
observations of the P V resonance doublet obtained with FUSE , ORFEUS/BEFS , and Copernicus. 

When P 4+ is the dominant ion in the wind (i.e., 0.5 < q , (P 4_t ') < 1), Mq( P 4+ ) approximates 
the mass-loss rate to within a factor of <2. Theory predicts that P 4+ is the dominant ion in 
the winds of 07-09.7 stars, though an empirical estimator suggests that the range from 04-07 
may be more appropriate. However, we find that the mass-loss rates obtained from P V wind 
profiles are systematically smaller than those obtained from fits to Ha emission profiles or radio 
free-free emission by median factors of ~130 (if P 4+ is dominant between 07 and 09.7) or ~20 
(if P 4+ is dominant between 04 and 07). These discordant measurements can be reconciled if 
the winds of 0 stars in the relevant temperature range are strongly clumped on small spatial 
scales. We use a simplified two-component model to investigate the volume filling factors of the 
denser regions. This clumping implies that mass-loss rates determined from “p 2 ” diagnostics 
have been systematically over-estimated by factors of 10 or more, at least for a subset of O stars. 
Reductions in the mass-loss rates of this size have important implications for the evolution of 
massive stars and quantitative estimates of the feedback that hot-star winds provide to their 
interstellar environments. 

Subject headings: stars: early-type - stars: mass loss - stars: winds, outflows 

1. Introduction 
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more, MD 21218. strong, dynamic influences on the evolution of 

massive stars and their interstellar environments. 
From the perspective of stellar evolution, contin- 
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uous mass-loss via these outflows alters both the 
path of the star through the H-R diagram and the 
rate at which it is traversed. The continual shed- 
ding of the outer layers of the atmosphere also 
causes systematic changes in the chemical compo- 
sition of the photosphere, particularly in the pres- 
ence of strong, rotationally-induced mixing. Ad- 
ditionally, the injection of chemically enriched ma- 
terial, momentum, and energy from the winds of 
hot, massive stars into their surroundings is an in- 
tegral part of the “stellar feedback” cycle, which 
mixes the local interstellar medium and ultimately 
drives the chemical evolution of a galaxy. 

In view of their fundamental astrophysical im- 
portance, it is clear that the properties of these 
outflows must be well determined. Nowadays, it is 
understood that hot-star winds are driven by the 
transfer of momentum from the stellar radiation 
held through scattering in metallic resonance lines. 
However, ab initio models still require numerous 
assumptions in order to provide testable predic- 
tions of wind parameters. Consequently, empirical 
confirmation of the theory remains essential. For 
this purpose, the three most useful diagnostics of 
mass loss are, in order of increasing sensitivity: 
(a) free-free continuum emission at radio wave- 
lengths; (b) Ho: line emission; and (c) ultraviolet 
(UV) resonance-line absorption. These diagnos- 
tics sample different parts of the wind, from the 
dense, near-star, rapidly accelerating region (Ha) 
to the very distant, rarefied, constant velocity re- 
gions (free-free radio emission) and essentially ev- 
erywhere between (UV resonance lines). Their 
physical mechanisms have different dependencies 
on the local density (emission oc p 2 ; absorption 
oc p), while the measurable properties associated 
with each suffer different degrees of contamination 
from plrotospheric or other processes, and require 
different ancillary knowledge of, e.g., the excita- 
tion or ionization conditions in the wind, the ve- 
locity structure of the wind, or the distance to the 
star. 

The need for additional information concerning 
the ionization structure of the wind has tradition- 
ally been a stumbling block for determinations of 
the mass- loss rates (M) from the wind profiles of 
UV resonance lines. The problem arises because 
the strength of these profiles is determined by the 
radial optical depth of the wind, r ra< i oc Mq^As, 
where Ae and qi are the abundance of element 



E and its ionization fraction for stage i. Thus, 
wind-profile modeling of a dominant ion (cp ~ 1) 
of known abundance is required to estimate M di- 
rectly. Unfortunately, the lines of dominant ions 
are usually saturated, particularly those for cosmi- 
cally abundant elements (e.g., C, N, O) in winds 
that are sufficiently dense to provide reliable de- 
terminations of M from “density squared” [here- 
after “p 2 ” and M(p 2 )) diagnostics. As a result, 
only UV wind lines from trace ionization species 
w ith Qi < 10 -3 can typically be used to measure 
M q for massive winds, and values of M estimated 
from these measurements are compromised by the 
lack of a priori information needed to estimate qi 
accurately; see Lamers et al. (1999) for additional 
discussion of this problem. Thus, even though UV 
wind lines are the most sensitive indicators of M in 
early-type stars, empirical mass-loss rates have re- 
lied primarily on measurements M(p 2 ) from free- 
free radio emission [M (radio)] and Ha line emis- 
sion [M(Ha)]; see, e.g., Lamers & Leitherer (1993) 
or Puls et al. (1996). 

Although mass-loss rates determined from ra- 
dio measurements are usually regarded as the most 
reliable, comparatively few OB-type stars have 
been detected owing to the inherent weakness of 
free-free emission and its “p 2 ” dependence, which 
requires stars to be relatively nearby. In con- 
trast, Ha emission, which is also a “p 2 ” processes, 
can be easily observed in distant objects. How- 
ever, because its strength depends on the details 
of the wind photosphere interface, it is more dif- 
ficult to model, although the current consensus is 
that mass-loss rates determined from Ha are quite 
accurate (Puls et al. 1996). 

Access to the far-ultraviolet region of the spec- 
trum provided by the Far Ultraviolet Spectroscopic 
Explorer (FUSE) affords a renewed opportunity 
to examine the use of resonance-line diagnostics 
for the determination of M, and to perform con- 
sistency checks on measurements of M(p 2 ). For 
this purpose, the Pv AA1118, 1128 doublet plays 
a pivotal role, primarily because P 4+ is expected 
to be the dominant ion the winds of some O stars. 
This expectation is based on naive energetic con- 
siderations, morphological trends observed in the 
spectra of O stars, and detailed modeling. From 
a diagnostic perspective, an equally important at- 
tribute of P 4+ is that, despite its high ionic abun- 
dance, the Pv AA1118, 1128 doublet is rarely sat- 
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urated because the cosmic abundance of P is so 
small. On the other hand, this also means that 
the wind contribution to the doublet is only mea- 
surable in stars with massive winds. Finally, in 
contrast to C, N, and O, P is not produced by H- 
burning; consequently its atmospheric abundance 
does not change substantially over a stellar life- 
time or as a function of the vigor of rotational 
mixing. Thus, star-to-star changes in its abun- 
dance should not contribute to the cosmic scatter 
in M determinations. 

These properties make the P V resonance dou- 
blet a powerful mass-loss diagnostic for a range 
of O-star temperature classes, though we expect 
the detection of wind profiles in the P V resonance 
doublet to be biased toward stars with more mas- 
sive winds. If we simply assume that g(P 4+ ) ap- 
proaches unity in stellar winds somewhere in the O 
star range, then values of Mq{ P 4+ ) should agree 
with both M (radio) and M (Ha) for at least some 
O stars. If concordance between these individu- 
ally reliable mass-loss indicators is not found, then 
their different formation mechanisms will help to 
constrain the reasons for the discrepancy. Massa 
et al. (2005) reported on a preliminary compari- 
son between these different measures of M , and 
the current paper expands and refines that work. 

In the present study, we have measured Mq(P 4+ ) 
for a sample of 40 well-studied Galactic O-type 
stars that have reliable radio or Ha mass-loss mea- 
surements available from the literature. This sam- 
ple is described in §2. New FUSE observations and 
wind-profile analysis of all P V resonance doublets 
are described in §3.1 and 4.3, respectively. The 
values of M q derived for P 4+ are compared with 
the values of M(p 2 ) in §5, and large discrepancies 
are found. Possible reasons for this discordance 
are discussed in §6, before the main conclusions 
from this study are summarized in §7. 

2. The Sample 

Table 1 lists the sample of Galactic O-type stars 
used in this investigation. The sample is defined 
by the availability of (a) reliably determined values 
of the fundamental stellar parameters; (b) reliable 
determinations of M from thermal radio emission 
or Ha emission profiles; and (c) far-UV spectra 
that include the P V resonance lines. For practical 
purposes, the first two criteria limit the selection of 



objects to those that have been well-studied over 
the past decade by the methods of “quantitative 
spectroscopy.” These objects serve as a fundamen- 
tal data set, and are the cornerstone for the cali- 
bration of the Wind-Momentum-Luminosity rela- 
tion for Galactic O stars (Puls et al. 1996; Repo- 
lust et al. 2004). The sample spans all subclasses 
of O-type spectra, but is biased toward higher lu- 
minosity classes. The luminosity bias occurs be- 
cause both “p 2 ” diagnostics and the P V wind fea- 
tures are preferentially detected in the densest out- 
flows, which are associated with more luminous 
objects. 

2.1. Stellar Parameters 

As described in §4.3, determinations of Mq( P 4+ ) 
from the P V resonance lines requires knowledge 
of the stellar radius, I?*; the projected rotational 
velocity, usini, which determines the width of the 
photospheric lines; and the terminal velocity of 
the stellar wind, Voq. These parameters and their 
origin are indicated in Table 1, along with the 
spectral type and the adopted values of T e g. 

To minimize possible systematic biases in the 
derived stellar parameters, the sample is subdi- 
vided into “primary” (28 objects) and “secondary” 
(12 objects) subsets according to how the parame- 
ters were determined. Most of the determinations 
for the primary data set come from the work of 
Repolust et al. (2004), which provides the largest 
sample of homogeneously determined stellar pa- 
rameters currently available for Galactic O stars. 
These are derived from fits made with the uni- 
fied (i.e. , photosphere + wind), non-LTE, line- 
blanketed model atmosphere program FASTWIND 
(Puls et al. 2005). The primary sample also con- 
tains objects analyzed by Markova et al. (2004) 
with a more simplified approach, which relies in 
part on calibrations of fundamental parameters 
derived from the work of Repolust et al. (2004). 
Thus, the values of the derived stellar parameters 
for the primary sample are internally consistent 
and incorporate recent revisions to the tempera- 
ture scale associated with O-type stars (Martins 
et al. 2002; Crowther et al. 2002; Bianchi & Gar- 
cia 2002; Herrero et al. 2002; Martins et al. 2005). 

Since similar fits with unified model atmo- 
spheres are not yet available for stars in the 
secondary sample, the values of T e g and R/Rq 
listed for them in Table 1 are based on the revised 
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spectral-type calibrations presented by Martins 
et al. (2005). The parameters given for the “ob- 
servational scale” (i.e. , Tables 4-6 of Martins et al. 
2005) were used to minimize possible systematic 
differences with the primary sample. Values for 
luminosity class II were estimated by linear in- 
terpolation between luminosity classes I and III; 
and the parameters associated with 09.5 super- 
giants were adopted for the three 09.7 supergiants 
in the secondary sample. Of course, it would be 
very useful to determine the stellar parameters of 
all the stars listed in Table 1 through rigorously 
uniform analysis. 

Table 1 also notes which objects in the sample 
are known or suspected binary stars, as designated 
in the Galactic O-Star Catalog (Mafz-Apellaniz 
et al. 2004). 

3. Far-Ultraviolet Spectra 

The far-UV spectra used to measure the Pv 
wind profiles of the targets in the sample were ac- 
quired by several observatories. These are indi- 
cated in Table 2, along with other salient details 
concerning the observations. Most of the spectra 
were obtained with FUSE, either as part of our 
Guest Investigator program E082, which was ded- 
icated to measurements of P V, or as a by-product 
of programs devoted to other topics. The process- 
ing of these data is described in §3.1. 

Owing to count-rate limits associated with the 
FUSE detectors, bright objects are difficult or im- 
possible to observe. Consequently, FUSE can- 
not observe nearby, unreddened O stars. This 
is unfortunate, since these are generally the ob- 
jects that have received the most attention from 
ground-based observers and for which radio mass- 
loss rates are available. For these bright stars, 
we have used archival far-UV spectra that were 
obtained by Copernicus or the Berkeley Extreme 
and Far-UV Spectrometer ( BEFS ). These spectra 
were retrieved from the Multi-Mission Archive at 
the Space Telescope Science Institute (MAST). 1 
The data sets for these instruments are indicated 
in Table 2, and were used without further manip- 

1 The archiving of non-HST data at MAST is supported by 
the NASA Office of Space Science via grant NAG5-7584 
and by other grants and contracts. STScI is operated by 
the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. 



ulation. 

3.1. FUSE Observations and Reduction 

The FUSE observatory consists of four aligned, 
prime-focus telescopes and Rowland-circle spec- 
trographs that feed two photon-counting detectors 
(Moos et al. 2000; Sahnow et al. 2000). These 
four channels provide redundant coverage of the 
range between ~ 905 — 1187 A. As indicated in 
Table 2, most objects were observed through the 
large 30" x 30"(LWRS) aperture. The lone excep- 
tion is HD 188209, which was observed through 
the narrow 1.25" x 20"(HIRS) aperture as part of 
a test of observing techniques for bright objects. 

The FUSE spectra were uniformly extracted 
and calibrated with CalFUSE version 2.4.2. Sub- 
sequent processing used the shifts determined by 
cross-correlating the positions of sharp interstel- 
lar features to align the spectra extracted from 
individual exposures; resampled the spectra onto 
a uniform wavelength grid with steps of 0.1 A; and 
normalized the profiles to the local stellar contin- 
uum. Since normalization of LiFIB spectra ac- 
quired through the LWRS aperture is complicated 
by an optical artifact known as “the worm” (Sah- 
now 2002), we used the redundant and cleaner 
coverage of the P V doublet provided by LiF2A 
spectra whenever possible. 

4. Mass-Loss Determinations 

Determinations of the mass-loss rates for O 
stars traditionally assume that the wind is a 
spherically symmetric, steady, homogeneous, and 
monotonically expanding outflow. All the mea- 
surements discussed in this section depend on 
these underlying assumptions of this “standard 
wind model.” 

4.1. Ha Measurements 

Most of the stars in the sample have M(Ha) 
measurements. The bulk of these were obtained 
from comparisons with FASTWIND calculations of 
the entire optical spectrum by Repolust et al. 
(2004), which were also used to derive other funda- 
mental stellar parameters. Since these computa- 
tions help to define the revised temperature scale, 
the values of M( Ha) are completely consistent 
with it. 
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An important subset of objects have lif(Ha) 
determined from fits to only the Ha line profile by 
Markova et al. (2004). They used the quick but 
accurate method developed by Puls et al. (1996). 
These measurements rely in part on calibrations 
of stellar parameters as a function of spectral type 
that were based on the complete modeling of Re- 
polust et al. (2004). Markova et al. (2004) show 
that their results are internally consistent with the 
values derived from the more complete analysis. 

Less reliable estimates of M (Ha) were used for 
four stars in the secondary sample: HD 45160, 
HD 46223, HD 164794, and HD 188001. For these 
stars, M was determined by Puls et al. (1996, 
their Table 11) from a re-analysis of the equiva- 
lent width measurements compiled by Laniers & 
Leitlierer (1993). These values are not of com- 
parable quality to the other determinations used 
here, since they (a) rely on equivalent widths 
rather than fits to the entire line profile; (b) are in 
some cases obtained from measurements of pho- 
tographic plates, which are inherently less precise 
than the CCD data used by Repolust et al. (2004) 
and Markova et al. (2004); and (c) were inter- 
preted on the basis of a temperature scale that 
is systematically hotter than the one used for the 
bulk of the stars analyzed here. Since three of 
these targets (HD 46150, HD 46223, HD 164794) 
only have upper limits for Mq( P 4+ ), the greater 
uncertainty inherent in their values of M(Ha) 
does not affect the comparison significantly. The 
fourth star, HD 188001, is the only program star 
in this group with detections of both M (Ha) and 
Mq( P 4+ ). Although the T e g of this star has been 
adjusted in Table 1, its value of M (Ha) has not 
been rescaled. 

4.2. Measurements of Radio Free- Free 
Emission 

Mass-loss rates determined from radio free-free 
emission are available for 40% of the stars in the 
sample. Most of these have been taken from the 
compilation of Lamers & Leitherer (1993), who 
provide a thorough discussion of this technique. 
These results are supplemented by the more re- 
cent detections of HD 14947 and HD 190429A by 
Scuderi et al. (1998). Following the procedure 
adopted by Lamers & Leitherer (1993), the uncer- 
tainties quoted by Scuderi et al. (1998) have been 
increased by ±0.10 dex to account approximately 



for uncertainties in the distance to the objects. 

However, all these determinations pre-date the 
revisions to the effective temperature scale for O 
stars. At first sight, M( radio) seems immune from 
these changes, since T e g enters the expression for 
M explicitly only very weakly as a logarithmic 
contribution to the Gaunt factor; see, e.g., equa- 
tion (5) of Lamers & Leitherer (1993). In fact, 
the resultant changes in My and stellar luminos- 
ity generally result in different estimates for the 
distance, d, to the various stars; see, e.g., Martins 
et al. (2005). Since M( radio) oc d 1,5 , a system- 
atic difference results between M (Ha) (which, ex- 
cept for HD 188001, are on the cooler temperature 
scale) and M(radio). 

To account for this dependence, we have sys- 
tematically adjusted the values of M(radio). For 
stars in the primary sample, we used the value of 
My quoted by Repolust et al. (2004) and Markova 
et al. (2004), together with photometry from the 
Galactic O-Star Catalog (Mafz-Apellaniz et al. 
2004), to estimate the distance modulus via the 
standard formula 

DM = V- M v - 3.1 E(B - V) (1) 

This leads to a revised estimate of the distance, 
d new , and a correction factor for the published 
value of M (radio) of f corr = (d new / doid) 1 ' 5 , where 
doid was listed by Lamers & Leitherer (1993) or 
Scuderi et al. (1998). The same approach was fol- 
lowed for stars in the secondary sample, though in 
this case internally consistent values of My were 
adopted from the spectral-type calibration of Mar- 
tins et al. (2005). Table 3 lists the correction fac- 
tors that were applied to make the measurements 
of M( radio) consistent with the measurements of 
Af(Ha). 

4.3. Pv Wind-Profile Fitting 

The strength of the absorption trough of a 
P Cygni wind profile at any normalized velocity 
w = v/voo (where Vao is the terminal velocity of 
the wind) depends on the radial optical depth of 
material along the line-of-siglrt at that velocity, 
T ra d(w), which is directly related to the column 
density. 

For the P v resonance lines (whose basic prop- 
erties are listed in Table 4), T ra d(w) for each pro- 
gram star was determined by using the “Sobolev 
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with Exact Integration” (SEI; Lamers et al. 1987) 
wind-profile fitting technique with the modifica- 
tions described by Massa et al. (2003). These 
modifications permit the relationship between op- 
tical depth and position in the absorption trough 
to be an arbitrary function, and also allow ap- 
proximately for the effects of interstellar absorp- 
tion, whenever required. The analysis proceeded 
in three steps. 

First, the parameters of the velocity law were 
defined. The velocity law, w(x), was parameter- 
ized in terms of the usual /3-law by 

w(x) = w 0 + (1 - l/x) 13 (2) 

where w = v/Vqo is the normalized velocity; x = 
r/i?* is the radial distance in units of the stellar 
radius; and where the radial gradient is controlled 
by the parameter /3, which is typically ~1. The 
adopted values of i>oo are listed in Table 1. These 
were determined from the shape of strongly satu- 
rated P Cygni profiles, typically the Civ AA1548, 
1550 doublet, and included an estimate of any ad- 
ditional velocity dispersion (often referred to as 
“turbulence”) in the wind, which is usually re- 
quired to reproduce the slope of the blue edge of 
the absorption trough. We found that adopting 
a value of 0.05 i>oo for this additional velocity dis- 
persion was sufficient to model the Pv profiles. 
Similarly, we adopted (3 = 1 for all initial fits to 
the P v doublet, and altered it by modest amounts 
only if significantly better fits resulted. The final 
values of (3 are listed in Table 5. 

Second, the Pv doublet in each spectrum was 
fitted with the adopted velocity law by adjusting 
the optical depth in the radial direction, T rad (w), 
in ~20 velocity bins distributed through the ab- 
sorption trough. As explained in detail by Massa 
et al. (2003) , determining the optical depth in ab- 
sorption in a given bin also determines the contri- 
bution to blue- and red-shifted emission, so that a 
satisfactory fit to the entire profile is constructed 
by stepping through the absorption trough from 
its blue edge to line center. Underlying photo- 
spheric profiles were approximated by Gaussians, 
but since their influence is confined to low veloc- 
ities (within iwsinf; see Table 1 for the adopted 
values), their inclusion has little effect on the de- 
termination of T ra d{w) for w > 0.2, except in cases 
where the wind profile is extremely weak. A bigger 
problem is the blend of Pv A1128 with a strong, 



excited transition of SilV in the spectra of later 
O-type stars. In these cases, greater weight was 
given to the Pv All 18 profile. Figures 1 and 2 
illustrate the quality of the fits for a selection of 
targets. 

Third, these determinations of r ra d(ic) were 
used to compute Mq(P 4+ ;w) for each star from 
the relation: 



Mq{ P 4+ ; 
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where f,j is the oscillator strength of the appro- 
priate component of the P V resonance line; Ap is 
the abundance of P relative to hydrogen by num- 
ber (Table 4); /r is the mean molecular weight of 
the plasma; and all other symbols have their usual 
meaning. In equation (3), we adopted p = 1.34 
(appropriate for a completely ionized plasma of 
solar abundance), together with the values of i?* 
and Vqo listed in Table 1. The measurements 
of Mqi(w) were averaged over normalized veloc- 
ities between 0.2 and 0.8 to obtain the values of 
Mq( P 4+ ) listed in Table 5. 

The errors associated with the derived T rad (w) 
depend on its magnitude. For T rad (w) < 0.05, the 
error can be as large as a factor of two. Similarly, 
for Trad(w) > 2.0 the error can also be large, while 
for Trad(vj) ~ 1.0, it is on the order of 10%. Con- 
sequently, we provide only upper limits for lines 
whose blue (stronger) component has an optical 
depth less than 0.05. None of the program stars 
has an optical depth in the red (weaker) compo- 
nent of the P v doublet greater than 2.0, so there is 
no need to define lower limits on T rad {w). Typical 
errors in parameters derived from the integrated 
optical depths, such as Mq(P 4+ ), are < 25%, and 
only quantities derived from the weakest lines have 
significantly larger errors. The value of f3 can also 
affect T rad (w). However, (3 is typically determined 
to within ±0.2, and uncertainties of this amount 
have very little effect on T rad {w) over the range 
0.2 < w < 0.8. 



5. Results 

5.1. Comparison of M(p 2 ) with Mq( P 4+ ) 

The published values of M( radio) [modified as 
described in §4.2] and M( Ha) are listed in Ta- 
ble 5. They are compared with the average values 
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of Mq( P 4+ ) in Figure 3 as a function of spectral 
type, measurement technique, and status as mem- 
bers of the primary or secondary sample. When 
more than one M(p 2 ) measurement was available, 
M (radio) was preferred, except in cases where the 
uncertainties associated the M( Ha) measurement 
were substantially smaller. Since only upper lim- 
its on both M{p 2 ) and Mq( P 4+ ) are available for 
HD 46150 and HD 217086, the values for these 
stars are not plotted in Fig. 3. The position of 
HD 149757 is not indicated either, since its very 
small value of Mq( P 4+ ) lies beyond the limits of 
the plot. 

Several things are apparent from Figure 3. 
First, Mq{ P 4+ ) is systematically smaller than 
M(p 2 ) by substantial amounts. Second, the size 
of the deviation does not depend on whether the 
fiducial M(p 2 ) was M (radio) or M( Ha). Third, 
any biases in the determinations of M for stars in 
the secondary sample are much smaller than the 
systematic deviations from the line indicating a 
1-1 correlation. Finally, a strong dependence on 
binarity is not evident. 

However, systematic trends in the distribution 
of deviations with spectral class are apparent. The 
sample can be divided into three parts: 

1. The mid O-type stars (04 - 07.5) with 
strong wind features (i.e., those with log 
Mq( P 4+ ) > —8) exhibit the smallest devia- 
tions from the values of M (Ha) or M (radio) . 

2. The earliest (02-03.5) and latest (08-09.7) 
O stars exhibit systematically larger devia- 
tions from the 1-1 correlation line. 

3. The largest deviations belong to a group 
of five mid-0 dwarfs and giants that have 
M(Ha) measurements, and only upper lim- 
its for Mq{ P 4+ ). These stars are (from 
left to right in Fig. 3: HD 42088 (06.5 
V), HD 46223 [04 V((f+))], HD 47839 [07 
V((f))] HD 217086 (07 Vn), and HD 203064 
[07.5 III:n((f))]. 

Figure 4 presents an alternate representation of 
this discrepancy in terms of an empirical estimate 
of the mean value of g(P 4+ ), 

q est =Mq(P 4+ )/M(p 2 ) . (4) 

In Fig. 4, q est is plotted as a function of T e g for dif- 
ferent luminosity classes. Although expected to be 



near unity over some range of T e g , it is never more 
than ~0.11 for supergiants, and becomes progres- 
sively smaller for less luminous stars (0.06 for gi- 
ants; 0.04 for dwarfs). Thus, M{p 2 ) is system- 
atically larger than Mq(Y i+ ), as also indicated 
in Fig. 3. Since the initial expectation was that 
Mq( P 4+ ) should be nearly equal to M{p 2 ) for at 
least some O-type stars, the size and systematic 
nature of the deviations requires explanation. 

5.2. For Which Stars is P 4+ Dominant? 

A crucial assumption implicit in our empha- 
sis on Pv is that g(P 4+ ) ss 1 for O stars span- 
ning some range of T e ff- Apart from naive ener- 
getic considerations, this assumption is supported 
by detailed model atmosphere calculations. In 
particular, smooth-wind models computed with 
FASTWIND indicate that this temperature range is 
~31 - 34 kK 2 , which corresponds to spectral types 
between 09.7 07.5 for supergiants and 09.5 - 

08.5 for dwarfs. 

However, the distribution of the empirical es- 
timates of g(P 4+ ) in Fig. 4 suggests that the ap- 
propriate range of T e g might be shifted to hotter 
temperatures. For more luminous objects, g e st ex- 
hibits a broad maximum between T e g ~ 34 - 40 kK 
(i.e., 07.5 - 04; mid-range O stars). A peak near 
~41 kK may also be indicated for dwarfs, though 
the behavior of g e st is not very clear for these ob- 
jects. 

A simple interpretation of the distribution of 
<7est with T e g is that it traces the rise and fall of 
P 4+ as the dominant ion in the wind. In this 
case, P 3+ would be the dominant ion for T e g < 
34 kK, while the balance shifts in favor of P 5+ for 
T eff > 40 kK. For both extremes, g(P 4+ ) <C 1, as 
implied by Fig. 4. However, this straightforward 
interpretation is complicated by the reliance of q es t 
on M{p 2 ). For example, biases in the determina- 
tion of M{p 2 ) will also bias g est . Evidently a more 
systematic modeling effort will be required to de- 
termine the temperature range for which P 4+ is 
the dominant ion. 

Irrespective of this uncertainty in the appropri- 
ate range of T e g, we conclude that there is a large, 
systematic discrepancy between Mq(V 4+ ) and the 
“p 2 ” determinations, M(Ha) and M(radio), for 

2 We are grateful to R.-P. Kudritzki and M. A. Urbaneja for 
providing us with a grid of unpublished FASTWIND models. 
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some O stars. 

1. If, as predicted by standard FASTWIND mod- 
els, P 4+ is dominant in the winds of 07.5 - 
09.7 stars, then Table 5 indicates that the 
median discrepancy for the subset of 15 lu- 
minous (i.e. , non-dwarf) stars with solid de- 
tections corresponds to M(p 2 )/Mq{ P 4+ ) = 
129, with minimum and maximum values of 
17 and 501 for HD 24912 and HD 209975, 
respectively. 

2. If, as suggested by a straightforward inter- 
pretation of the behavior of q e s t in Fig. 4, 
P 4+ is dominant for stars with spectral 
classes between 04- 07.5, the median dis- 
crepancy for the subset of 13 luminous 
stars with solid detections corresponds to 
M(p 2 )/Mq( P 4+ ) = 20, with minimum and 
maximum values of 9 and 245 for HD 66811 
and HD 15558, respectively. 

5.3. Stars with weak winds 

The extremely large discrepancies for the mid 
O-type stars with weak winds (i.e., the class III 
- V stars) are an interesting special case, which 
might result from systematic measurement errors 
in M{ Ha). None of these targets have obvi- 
ous wind profiles in P V, so only upper limits for 
Mq{ P 4+ ) can be estimated. Similarly, none of 
these stars are detected in the radio, and none ex- 
hibit bona fide Ha emission profiles. As a result, 
these determinations of M (Ha) depend sensitively 
on estimates of the degree to which the under- 
lying photospheric profiles are partially filled by 
wind emission, which in turn relies critically on the 
accuracy of the photospheric models. Such wind 
contamination is clearly evident in the Ha pro- 
files of HD 203064 (Repolust et al. 2004), though 
emission from a circumstellar disk is not precluded 
for this rapid rotator. HD 217086 is also a rapid 
rotator, though in this case partial filling by emis- 
sion is not evident, and only an upper limit on M 
was determined (Repolust et al. 2004). Emission 
is less evident in the Ha profiles of HD 47839, 
which might instead be weakly contaminated by 
the spectrum of its late-0 companion (Gies et al. 
1993). 

In any case, if the determinations of M(Ha) 
for the less-luminous, mid O-type stars are taken 



at face value, they yield values that are inconsis- 
tent with the non-detection of wind profiles in the 
resonance lines of P V, at least for the assumed 
solar abundance of P. Although the magnitude of 
the discrepancy is indeterminate for these objects, 
the sense is the same as for the rest of the mid- 
O stars. The largest value of Mq{ P 4+ ) for any 
main sequence O star in our sample is 7.2 x 10~ 8 
M©/yr, and most are less than this by an order 
of magnitude or more. Thus, if g(P 4+ ) ever ap- 
proaches unity for main sequence O stars, then 
their mass-loss rates are much smaller than theo- 
retical expectations, which are > ICR 7 Mg/yr for 
these stars (Vink et al. 2000). 

5.4. Comparison of M(radio) with M(Ha) 

Although not a primary motivation of the 
present study, the compilation of measurements of 
M data in Table 5 provides an opportunity to re- 
examine the agreement between measurements of 
M (radio) and M (Ha) . Detections are available for 
both diagnostics for a subset of 8 objects. A com- 
parison of the two measurements for these stars in- 
dicates a systematic difference between them, with 
the mean value of Af(Ha) / M{ radio) « 2.4 ±0.7. 
This is contrary to previous studies (e.g., Larners 
& Leitherer 1993; Puls et al. 1996), which found 
agreement between the two measures. Blonnne 
et al. (2003) noted a similar - though smaller 
discrepancy between M (radio) and M (Ha) deter- 
minations for HD 66811 [04 I(n)f], 

6. Discussion 

It has been known for some time (see, e.g., 
Hamann 1980) that self-consistent wind models 
produce too much absorption in the P V resonance 
lines of individual O-types stars. Although the 
precise range of spectral types where g(P 4+ ) is 
dominant is poorly defined at present, the large 
sample considered here indicates that there is a 
significant, systematic discrepancy between mass- 
loss determinations for all objects with spectral 
types in either of the two possible ranges discussed 
in §5.2. Thus, the “Pv problem” is not limited 
to a few, possibly peculiar objects. Instead, the 
discordance between M diagnostics indicates that 
one or more biases must be influencing the anal- 
yses in such a way that either Mq( P 4+ ) is sys- 
tematically under- estimated, or that M( radio) and 




M(Ha) are systematically over-estimated. A com- 
bination of biases might also be affecting these di- 
agnostics by different amounts. In the following 
sections, we examine the biases that might affect 
mass-loss rates determined by different methods 
within the context of the standard wind model 
(i.e., a smooth, steady, and spherically symmet- 
ric wind). 

6.1. Biases in Measurements of M(Ha) 

A measurement of M{ Ha) is complicated, be- 
cause it requires (a) accurate knowledge of the 
distribution of wind density in the near-star envi- 
ronment; (b) non-LTE calculations to determine 
the excitation equilibrium of H; and (c) an accu- 
rate model for the underlying plrotospheric fea- 
ture, which must also incorporate the blend with 
He II A6560. The latter two requirements are both 
sensitive to temperature, with hotter tempera- 
tures usually leading to larger values of M (Puls 
et al. 1996). Consequently, values of M{ Ha) may 
be systematically over-estimated with old, hotter 
temperature scale. However, since only one of the 
older M( Ha) measurements was used in Figs. 3 
and 4 (see §4.1), this is unlikely to be an impor- 
tant bias. If anything, the reduction in the tem- 
perature scale should move values of M (Ha) closer 
to Mq{ P 4 +). 

A more subtle problem occurs for stars of lumi- 
nosity class V— III, which typically exhibit a par- 
tially filled absorption profile rather than emis- 
sion above the local continuum. For these stars, 
estimates of M(Ha) depend very sensitively on 
how well the underlying photospheric profile is re- 
produced by the adopted model atmosphere. It 
is possible that small, systematic uncertainties in 
the strength of the photospheric Ha profiles con- 
tribute to the large discrepancy between Mq { P 4+ ) 
and M (Ha) exhibited by the group of 5 mid-range 
O dwarfs and giants seen in Fig. 3. 

Another potential source of bias in determina- 
tions of M (Ha) is the high frequency of variability 
(Ebbets 1982; Kaper et al. 1998; Markova et al. 
2005). Many of these variations are characterized 
by changes in the equivalent width of the emission, 
which implies that either the amount of material 
or its density distribution changes. Consequently, 
it is difficult to evaluate the extent to which a mea- 
surement based on a single profile at a snapshot 
in time represents the global, tinre-averaged mass 



flux. Wind-wind interactions in short-period bi- 
nary systems also produce complicated variations 
in Ha; see, e.g., the series of papers bounded by 
Gies & Wiggs (1991) and Thaller et al. (2001). 
Although the effects of colliding winds are not in- 
cluded in the analyses used to compile Table 5, 
there is no particular evidence that that binaries 
depart from the general trends evident in Fig. 3. 

6.2. Biases in Measurements of Af (radio) 

The M (radio) are usually assumed to be the 
most reliable, since the physical mechanism re- 
sponsible for the radio emission (free-free emis- 
sion) is well understood; ionization corrections are 
generally insignificant (Laniers & Leitherer 1993); 
and contamination from the stellar photosphere 
is completely negligible. Furthermore, the radio 
photosphere is typically sufficiently far from the 
star that the details of the velocity law are not 
important: only i>oo is required, which can be de- 
termined accurately from UV resonance lines. In 
the context of the standard wind model, the only 
other input to M(radio) is the distance. Although 
imperfect knowledge of the distances to Galactic O 
stars introduces scatter into the derived M (radio), 
the uncertainties should be random and are un- 
likely to result in a systematic bias. 

However, it has been recognized for many years 
that the M (radio) can be systematically over- 
estimated if a fraction of the observed radio flux is 
due to nonthermal emission. Approximately 25% 
of O-type stars exhibit nonthermal radio emis- 
sion, which can be recognized either through dra- 
matic variability between epochs or from the fre- 
quency dependence of their radio spectra (Bieging 
et al. 1989). Binarity may play a role in gener- 
ating nonthermal electrons through wind-wind in- 
teractions; see, e.g., the recent discussion of ra- 
dio observations of HD 93129A by Benaglia & 
Koribalski (2004). However, other mechanisms 
may contribute to the nonthermal radio emission 
from the winds of single stars; see, e.g., Van Loo 
et al. (2004). Of the 16 stars in our sample with 
M (radio) measurements, only 3 are considered 
“definite” thermal emitters (HD 66811, HD 37742, 
and HD 152408) by Lamers & Leitherer (1993), 
and 2 more are considered “probable” thermal 
emitters (HD 30614 and HD 151804). The re- 
maining 11 stars have measurements at a single 
frequency, which is insufficient to establish that 
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the emission is thermal. Five of these objects 
are known or suspected binaries, although these 
stars do not appear to have systematically larger 
M (radio) values compared to similar, single stars. 
Most of the stars with M(radio) measurements 
also have M( Ha) measurements that are larger 
still; see §5.4. Thus, it seems that any biases com- 
mon to “p 2 ” diagnostics affect M{ Ha) more than 
M (radio). 

6.3. Biases in Measurements of Mq( P 4+ ) 

Equation (3) shows that Mq{ P 4+ ) will be sys- 
tematically under-estimated if either the oscillator 
strength, f t j , or the P abundance, Ap, is over- 
estimated. Although random errors in the other 
parameters required to evaluate equation (3) will 
introduce uncertainties, the resulting effects are 
unlikely to be systematic. 

Since P 4+ is a lithium-like ion, it is quite un- 
likely that its oscillator strength is uncertain by 
a substantial amount. Morton (2003) does not 
indicate that there is any controversy related to 
the value of f,j . Consequently, there is no rea- 
son to suspect that the oscillator strength biases 
measurements of Mq{ P 4+ ). 

The Galactic P abundance is more problem- 
atic. Some authors (e.g., Pauldrach et al. 1994, 
2001) have adopted subsolar values (e.g., P/P© of 
0.50-0.67 and 0.05 for HD 66811 and HD 30614, 
respectively) in order to achieve good fits for wind 
profiles in the Pv resonance doublet. However, 
there is little evidence to support a systematically 
reduced abundance of P for Galactic O-type stars. 
The solar P abundance by number relative to H is 
A P = 12.00 + log(AT P /A^) = 5.45±0.06 (Biemont 
et al. 1994)), which is ~30% smaller than the me- 
teoric abundance ( Ap = 5.56; Anders & Grevesse 
1989). Dufton et al. (1986) measured Ap = 5.59 
in the Galactic interstellar medium from a survey 
of the Pii AA1153, 1302 resonance doublet along 
51 sight lines and found that it was only mod- 
estly depleted (factors of 3 or less) in cold clouds. 
More recently, Lebouteiller et al. (2005) have also 
confirmed that the interstellar abundance of P is 
solar. Similarly, analysis of a FUSE spectrum of 
HD 207538 (B0 V, though a candidate for a chem- 
ically peculiar star) by Catanzaro et al. (2003) 
shows that the plrotospheric lines of P IV are well 
reproduced with solar abundances. Thus, there is 
no particular evidence that the abundance of P is 



systematically subsolar in the material from which 
O stars form, or in their cooler, B-type relatives. 
It is in principle possible to search for anomalous 
P abundances more explicitly by determining the 
plrotospheric abundance of P V from O-type stars 
with weak winds (e.g., HD 47839; see Fig. 2). 

A more speculative explanation for the system- 
atically small values of Mq( P 4+ ) is that g(P 4+ ) 
peaks for O stars at a value that is substantially 
less than 1. Of course, this behavior is not pre- 
dicted by standard wind models, and would re- 
quire that some unusual circumstances dominate 
the ionization balance of P. A candidate mecha- 
nism noted by Pauldrach et al. (1994) is the ap- 
proximate coincidence of the ground-state ioniza- 
tion threshold of P V (65.023 eV; 191 A) with that 
of Hen (54.416 eV; 228 A). As a result, Pauldrach 
et al. (1994) speculated that g(P 4+ ) is driven by 
the behavior of He + . To the best of our knowl- 
edge, this idea has not been developed further. 
Another candidate mechanism is the production of 
soft X-rays, which is already beyond the strict con- 
straints of the “standard model,” but which could 
alter the ionization balance in the wind. However, 
there is no evidence that ions with observable wind 
lines and ionization states that lie above and be- 
low P 4 + - e.g., S 5+ and P !+ - are strongly over- 
populated; see the Appendix A for a plausibility 
argument in the case of HD 66811. 

We conclude that, in the context of the stan- 
dard wind model, there is no compelling reason 
why measurements of M g(P 4+ ) should be system- 
atically under-estimated by factors of > 10. 

6.4. Relaxing the Standard Model 

The preceding discussion does not clearly iden- 
tify any processes that could account for the sys- 
tematic discrepancy of M measurements made 
from what should otherwise be reliable diagnos- 
tics. Instead, it would seem that the discordance 
is an artifact of one or more of the assumptions 
inherent to the “standard wind model.” 

The weakest assumption is that the wind has 
a smooth density distribution. The rampant vari- 
ability of UV and Ha wind profiles (Prinja & 
Howarth 1986; Kaper et al. 1996, 1998; de Jong 
et al. 2001; Markova et al. 2005); extended “black 
troughs” and variable blue wings in saturated ab- 
sorption troughs of P Cygni profiles (Lucy 1982; 
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Puls et al. 1993); stochastically variable sub- 
structure in the emission-line profiles of HD 66811 
(Eversberg et al. 1998); and the detection of X- 
rays distributed through the wind (Harnden et al. 
1979; Seward et al. 1979; Chlebowski et al. 1989; 
Cassinclli et al. 2001) all denote the presence of 
inhomogeneous structures on a variety of spa- 
tial scales. Theoretical calculations (Owocki et al. 
1988; Feldmeier et al. 1997) also indicate that line- 
driven stellar winds are subject to strong instabil- 
ities, which should redistribute wind material into 
dense structures. 

6. 4-1- The effect of clumping on “p 2 ” diagnostics 

The presence of clumping substantially alters 
the interpretation of “p 2 ” diagnostics like free-free 
continuum emission and Ha line emission. Since 
the emission results from the interaction of two 
particles, it will be produced more strongly from 
denser regions. Consequently, if a clumped wind 
is interpreted in terms of the “standard model” 
(which is smooth and homogeneous), “p 2 ” di- 
agnostics will necessarily over-estimate the true 
mass flux, since the excess emission produced in 
the inhomogeneities will be incorrectly interpreted 
as arising from a smooth but denser medium. 
Smooth but aspherical redistributions of wind ma- 
terial such as those caused by rapid rotation (see, 
e.g., Petrenz & Puls 1996) also lead to over- 
estimates of the mass-loss rate by “p 2 ” diagnos- 
tics, again because the standard model provides 
an incorrect density distribution. 

Of course, the sensitivity of “p 2 ” diagnostics 
to clumping has been recognized for a long time. 
It had been previously discounted on the basis 
of the perceived agreement between M(Ha) and 
M (radio), which was taken to indicate that either 
the wind is clumped by the same amounts over 
huge distances (which was viewed as unlikely) or it 
is not significantly clumped anywhere (Laniers & 
Leitherer 1993). The reassessment of these mass- 
loss rates, which is partially driven by the lowering 
of the T e s scale for O-type stars, has challenged 
the observational basis for this conclusion (Repo- 
lust et al. 2004). 

Clumping is frequently described in terms 
a two-component model that consists of dense 
clumps characterized by density pn and a rar- 
efied inter-clump gas with density pl , which re- 
distributes the material of a smooth flow while 



preserving its mean density, p s : 

Ps = f pH + (1 - /) PL = PH {f + [1 - /] x) (5) 



where / is the volume filling factor of the denser 
component, 0 < / < 1; x is the density con- 
trast, x = pl/ph'i and where all quantities are 
understood to be functions of position. Although 
initial applications of this formulation empha- 
sized continuum processes (see, e.g., Abbott et al. 
1981), it has also been used characterize the ef- 
fect of clumping on line processes (Hillier & Miller 
1998; Crowther et al. 2002; Bouret et al. 2005). 
Since the physical mechanism responsible for re- 
distributing wind material into dense clumps is 
not known, most implementations minimize the 
number of free parameters by further assuming 
that the dense clumps are separated by vacuum, 
so that x = 0 and p s = f pH- To put the present 
discussion on the same footing, we also follow this 
approach. 

As shown by Abbott et al. (1981), the two- 
component model for clumping predicts that the 
M determined from measurements of “p 2 ” diag- 
nostics, M(p 2 ), will be over-estimated by a factor 
of 



( {f + x(l-f)] 2 \ h 

\ / + (1 - f)x 2 J 



( 6 ) 



In the extreme case of dense clumps separate by 
vacuum, x = 0 and M is over-estimated in a 
smooth- wind model by i.e. , M(p 2 ) c = ffj 

M(p 2 ) s , where the subscripts “c” and “s” de- 
note “clumped” and “smooth,” respectively. In 
effect, this simplified model for clumping predicts 
a degeneracy between f and M(p 2 ); see Hillier & 
Miller (1999). 



6-4-2. The effect of clumping on UV resonance 
lines 

In contrast to “p 2 ” diagnostics, determinations 
of Mq from the wind profiles of UV resonance lines 
are quite insensitive to clumping. The analysis of 
these profiles is essentially a determination of the 
optical depth (or, equivalently, the column den- 
sity) of all the material associated with a specific 
ion along the line of sight that causes the observed 
P Cygni absorption trough. Since optical depth 
(or column density) is an integral quantity, these 
measurements are not sensitive to the distribution 
of material along the line of sight. 
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Consequently, as long as the clumps remain op- 
tically thin on the spatial scales relevant to line 
transfer, no material will be hidden and measure- 
ments will reliably account for all of it irrespective 
of its distribution. Thus, from equation (3), ( Mq) s 
= ( Mq ) c . The situation becomes more compli- 
cated if the clumps become optically thick, since 
their shape and distribution must then be specified 
in order to determine the degree to which the wind 
is “porous.” Owocki et al. (2004) provide a useful 
formalism to describe porosity in this context, and 
discuss how it can affect the predicted mass-loss 
rates. Massa et al. (2003) describe how porosity 
can affect the formation of P Cygni lines and, in 
extreme cases, produce an apparently unsaturated 
profile for a line that would be extremely saturated 
if the wind material were distributed smoothly. 

On larger spatial scales, inhomogeneities will 
be directly observable as significant departures 
from the expected shape of the P Cygni absorp- 
tion trough. The discrete absorption components 
(DACs) that are nearly ubiquitous in wind profiles 
of O-type stars are notable examples of such fea- 
tures; see, e.g., Prinja & Howarth (1986); Howarth 
& Prinja (1989); Kaper et al. (1996, 1999). The 
SEI-fitting procedure described in §4.3 will char- 
acterize the optical depth of large-scale structures 
like DACs reliably. 

The detailed model-atmosphere analysis of 
HD 190429A (04 If+) by Bouret et al. (2005) 
provides a test of the robustness of SEI fits to 
small-scale clumping. Our SEI fit to the P V pro- 
files in this star is excellent (Fig. 1), and conse- 
quently Mq{ P 4+ ) is very well determined. Simi- 
larly, Bouret et al. (2005) achieved an excellent fit, 
but only for a “clumped” model with a P abun- 
dance reduced to P/P© = 0.5; see their Fig. 2. 
Using this value of P/P©, their M (from their 
Table 4), and their mean g(P 4+ ) = 0.5 over the 
velocity range 100-2000 kins -1 (estimated from 
their Fig. 5), we find that the clumped CMFGEN 
model predicts logM< 7 (P 4+ ) = —6.08. For com- 
parison, scaling the results of our SEI profile fitting 
by their reduced P abundance gives logMg(P 4+ ) 
= —6.06. This remarkable agreement confirms 
that both techniques determine the same optical 
depth in the line when clumping is incorporated 
in the models and all other factors are equal; i.e. , 
clumping does not bias the determination of Mq 
from wind-profile fits to UV resonance lines. 



However, even though the robustness of wind- 
profile fitting against clumping ensures that Mq 
will be determined reliably, we expect that q will 
be affected by the presence of significant clumps 
in the wind. Two-body interactions occur more 
frequently in denser clumps, so recombination to 
lower states is favored. In general, 
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The ratio Q = g s (P 4+ )/g c (P 4+ ) can be greater 
or less than 1, depending on whether q s is domi- 
nant ( q c decreases with respect to q s ; Q > 1 ) or 
one stage below dominant ( q c increases with re- 
spect to q s ; Q < 1). For example, the models of 
HD 190429A by Bouret et al. (2005, (their Fig. 5) 
show that in the high-velocity region of the wind, 
q s « 0.1 and q c « 0.6; i.e., Q « 0.16. Q can in 
principle be exactly 1 if clumping causes the gain 
in g(P 4+ ) from P 5+ to balance the loss to P 3+ . 
Unfortunately, it is not possible to determine the 
relevant populations of these ions in a purely em- 
pirical fashion, though we show in Appendix A 
that some constraints can be inferred for the spe- 
cial case of HD 66811 (/ Puppis). 



6.5. Reconciling Mass-Loss Estimates 

The different behavior of “p 2 ” diagnostics and 
the resonance lines of dominant ions in the pres- 
ence of small-scale clumping implies that concor- 
dance between these measurements can in princi- 
ple be achieved by allowing for inhomogeneities 
in the wind. In particular, measurements of 
M(Pv) can be used to break the degeneracy be- 
tween M(p 2 ) and / for stars where g(P 4+ ) ~ 
l. 3 In terms of the simplified model for clump- 
ing discussed in §6.4.1, concordance between these 
two estimates can be achieved when M(Pv) c = 
M(p 2 ) c ; i.e., when the filling factor is 



/ 



Q 2 



( m(pv) s \ 

V M(P 2 )s ) 



( 8 ) 



Unfortunately, Q is a function of both the degree 
of clumping and T e ff- Although this information is 
well defined in terms of the parameters of specific 



'^Bouret et al. (2005) noted that the radiatively pumped 
O IV AA1338, 1343 and O V A1371 wind lines also help to 
break the M (p 2 )—f degeneracy. 



12 




model atmosphere computations, its behavior is 
difficult to predict a priori. Thus, for illustrative 
purposes, we arbitrarily assign Q values between 
0.1 and 10. We expect this range to cover most 
cases where a P V wind profile is observed, irre- 
spective of whether smooth wind models predict 
P 4+ to be dominant (q s « 1; Q > 1 ) or a level 
below dominant (q s <C 1; Q < 1). 

Table 6 shows the filling factors implied by 
these assumptions for the two possible ranges of 
T e ff discussed in §5.2. If smooth wind models cor- 
rectly predict the T e s range where < 7 S (P 4+ ) « 1, 
then Q > 1; but if the hotter range indicated by 
q es t is appropriate, then Q < 1. In either case, Ta- 
ble 6 shows that substantial degrees of clumping 
are required to achieve concordance, irrespective 
of the adopted value of Q. 

The filling factors associated with the simpli- 
fied model of clumping introduced in §6.4.1 have 
also been determined from comprehensive spec- 
troscopic analysis of individual objects. To date, 
these analyses have been undertaken with CMFGEN 
(Hillier & Miller 1998), which implements velocity- 
dependent filling factors that decrease exponen- 
tially from 1 (i.e., smooth) near the photosphere to 
a “terminal” filling factor, /<*,, at great distances 
from the star. Since the exponential decrease in / 
is very rapid, the values of foo from these analyses 
should correspond to the mean values of / implied 
by the mass-loss discrepancy. The values of 
derived from detailed spectroscopic fits to O stars 
in the Small Magellanic Cloud are typically ~0.1 
(Crowther et al. 2002; Hillier et al. 2003; Bouret 
et al. 2003; Evans et al. 2004). Bouret et al. (2005) 
found that clumped models with values of foo = 
0.04 and 0.02 were required to fit the UV spectra 
of two Galactic O stars, HD 190429A (04 If+) 
and HD 96715 [04 V((f))], respectively. Although 
these values of / are not as extreme as the illustra- 
tive values provided in Table 6, we conclude that 
qualitatively, both detailed spectroscopic model- 
ing and a straightforward analysis of the mass- 
loss discrepancy indicate that the winds of at least 
some O stars are strongly clumped. 

Of course, it remains to be seen whether such 
small filling factors are physically realistic. The 
application of the simplified model of clumping 
described in §6.4.1 to line transfer has not been 
justified rigorously. Although it likely captures 
some of the effects associated with density en- 



hancements, it is fundamentally inappropriate for 
line processes because it does not specify a length 
scale. As a result, the clumps are assumed to re- 
main optically thin irrespective of their density en- 
hancement, which is unlikely to be appropriate for 
very small filling factors (which imply very large 
values of pn/ps)- In addition, the usual approx- 
imation for the density contrast (i.e., x = 0) is 
aphysical (“Nature abhors a vacuum” 4 , etc.), and 
leads to extreme values of /. Finally, observa- 
tions of stellar-wind profiles indicate that large- 
scale structures are generally present, and these 
are likely to contribute their own perturbations 
to the assumed density structure in addition to 
small-scale clumping. For all these reasons, the 
filling factors derived from the simplest version of 
the “two-component” model are likely to be too 
small. 

7. Conclusions 

Comparison of Mq{ P 4+ ) values determined 
from wind-profile fitting with mass-loss rates de- 
termined from radio and Ha emission reveals a 
large, systematic discrepancy between diagnostics 
that should be reliable mass-loss estimators. The 
magnitude of this discrepancy depends crucially 
on the range of T e g for which P 4+ is the dom- 
inant ion in the wind. If, on the one hand, the 
predictions of the ionization balance of P based 
on smooth-wind models are adopted, then P 4+ is 
dominant between 07.5 and 09.7. In this case, the 
median discrepancy is enormous: M(p 2 )/M( Pv) 
« 130. If, on the other hand, the “p 2 ” diagnostics 
are taken at face value, then g est (P 4+ ) is observed 
to peak for temperature classes between 04 and 
07 at values of ~11% for supergiants, 6% for gi- 
ants, and ~4% for dwarfs. The fact that these 
maximal values are not unity underscores the dis- 
cordant estimates of M . The median discrepancy 
over this subset of our sample is more modest: 
M(p 2 )/M( Pv) ss 20. Massa et al. (2003) found 
similar behavior for q e st(P 4+ ) hr a sample of O 
stars in the Large Magellanic Cloud. 

Thus, irrespective of the range of T e g for which 
g(p 4 +) ~ i, we conclude that there is a signifi- 
cant discrepancy in M estimates for some Galactic 
O stars. We interpret the disagreement between 

4 Attributed to Frangois Rabelais (c. 1494-1553), a French 
monk, satirist, and physician. 
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these diagnostics as a signature of the presence of 
significant, small-scale clumping in the winds of 
some O-type stars (and, by likely extension, all O 
stars, though perhaps to substantially different de- 
grees). If the winds are clumped, then the values 
of M measured from “p 2 ” diagnostics are over- 
estimated, while the values of Mq measured from 
a dominant ion like P 4+ will be modestly affected, 
if they are affected at all. On the basis of a simple 
model for clumping, we find that consistent esti- 
mates of the mass-loss rate can be recovered if the 
dense material is characterized by volume filling 
factors that are very small. 

Since previous determinations and calibrations 
of M have relied on “p 2 ” diagnostics, we conclude 
that mass-loss rates have been over-estimated for 
at least some O-type stars. Thus, the primary con- 
sequence of the present study is that estimates of 
the mass-loss rate must be decreased substantially, 
typically by factors of ~ 100 [if <?(P 4+ ) ~ 1 between 
07.5 and 09.7] or ~20 [if g(P 4+ ) ~ 1 between 
04 and 07]. Even larger decreases are implied 
for mid-range dwarfs. Interestingly, independent 
studies of the X-ray emission from O stars - in par- 
ticular HD 36486A (6 Orionis A; 09.5 II; Miller 
et al. 2002) and £ Pup (Kramer et al. 2003) - have 
found that the cool wind contributes significantly 
less attenuation than expected for the adopted 
mass-loss rates, which also suggests that M has 
been over-estimated. Revisions to M of the mag- 
nitude suggested by the present study will have 
wide-ranging consequences for both the evolution 
of massive stars and the feedback they provide to 
their interstellar environments. 

Our conclusions concerning the importance of 
clumping are similar to those reached in studies of 
the winds of Magellanic O stars (Crowther et al. 
2002; Massa et al. 2003; Hillier et al. 2003; Bouret 
et al. 2003; Evans et al. 2004); a survey of lumi- 
nous Galactic B stars (Prinja et al. 2005); and the 
detailed analysis of two Galactic 04 stars (Bouret 
et al. 2005). All these studies highlight the util- 
ity of unsaturated resonance lines of dominant or 
near-dominant ions as mass-loss indicators. The 
P V resonance lines are the crucial diagnostics for 
O-type stars; Evans et al. (2004) have suggested 
that the Siv AA1073, 1074 doublet plays an anal- 
ogous role for spectroscopic analysis of late-0 and 
early B-type supergiants. Access to these lines 
must now be considered essential for detailed spec- 



troscopic studies of early-type stars. 

At an even more fundamental level, this work 
indicates that the “standard model” for hot-star 
winds needs to be amended, since predictions 
based on it appear to be inadequate for many 
purposes. In particular, clumping needs to be in- 
corporated realistically as a general feature. In 
order to make this feasible, the mechanisms re- 
sponsible for the formation of clumps must to be 
understood, as must the evolution of these inho- 
mogeneities under the ambient radiative and hy- 
drodynamic forces. The implication is that knowl- 
edge of the time-dependent behavior of line-driven 
stellar winds on large- and small-spatial scales is 
required to understand even the basic properties 
of these outflows, and in particular to interpret 
their observational signatures correctly. 

We wish to thank the referee, Rolf Kudritzki, 
for his detailed comments on this paper. Many 
improvements resulted from his insightful criti- 
cisms. We are likewise grateful to Joachim Puls for 
fruitful discussion of many of the issues explored 
here. This work is based in part on observations 
made with the NASA-CNES-CSA Far Ultraviolet 
Spectroscopic Explorer, which is operated by the 
Johns Hopkins University for NASA under con- 
tract NAS5-32985. Financial support from NASA 
Research Grants grants associated with FUSE 
Guest Investigator program E082 (NNG04EH17P 
to SGT, Inc. and NNG04GM49G to JHU) is also 
gratefully acknowledged. 

Facilities: FUSE, Copernicus, Orfeus(BEFS). 
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A. Empirical Constraints on the Ionization of P in the Wind of HD 66811 

Ideally, the ionization state of P in a stellar wind could be probed directly through the ratios 
g,(p5+)/^(p 4 +) anf ] q(P 3 +)/g(P 4 +). Unfortunately, the resonance lines of P 5+ lie near 90 A and are 
inaccessible. The resonance line of P 3+ falls at 951 A, which is accessible to FUSE but almost always 
compromised by extinction or blending; see, e.g., the far-UV atlases of Pellcrin et al. (2002) and Walborn 
et al. (2002). 

Instead, S 3+ and S 5+ can be used as surrogates for P 3+ and P 5+ to constrain the ionization state of 
P. Although these matches are not perfect, they are useful indicators of the prevailing ionization conditions 
in the winds of mid-0 stars, because (a) the ionization ranges of the surrogates correspond approximately, 
and in any case bracket, the range of P 4+ (see, e.g., Massa et al. 2003, their Fig. 1); and (b) like P, S is a 
non-CNO element whose abundance does not change over the life of a given star, and is unlikely to exhibit 
gross abundance differences between stars. Owing to its modest reddening, HD 66811 [<( Puppis; 04 I(n)f] 
is one of the few Galactic O stars for which the S VI resonance doublet can be measured. Consequently, the 
following plausibility argument is based on the Copernicus spectrum of £ Pup (Morton & Underhill 1977). 
The salient characteristics of this spectrum are that the S iv wind lines are weaker than those of P V, while 
the wind lines of S VI are strong but not saturated. 

To connect the ionization properties of P and S, consider that for two wind lines of the same star (where 
the iv(x) and M are identical for all lines), equation (3) gives 



q W 
g( 2 ) 



r (1) 

rad 

r (2) 

rad 



A f/ifU 2 ' 

A 0 Jij a e 



(Al) 



Since the products of XofijAp; are known for each resonance line and the values of T ra d can be determined 
by wind-profile fitting, the relative strengths of wind lines can be used to infer the relative ion fractions. For 
solar abundances, the ratios of Aq A e for Siv/Pv and Svi/Pv are approximately 6 and 59, respectively. 
Then, by using the observed S ions as surrogates for the unobserved P ions: 



1. The observation that the Siv doublet is weaker than the Pv doublet in the spectrum of HD 66811 
implies that T ra< i(Slv)/T ra( i(P v) < 1, and hence q( P 4+ ) > 6q(S 3+ ) « 6g(P 3+ ). 

2. The fact that the Svi doublet is not saturated implies that r ra d(Svi) < 5, while T ra d(Pv) ~ 1. Hence 
g(P 4+ ) > 12<AS 5+ ) » 12g(P 5 +). 

Thus, in the wind of HD 66811, we infer that q ( P 4+ ) is greater than the ion fractions of either of its adjacent 
states. While this does not in itself prove that g(P 4+ ) «1, it shows that it is not possible to hide significant 
amounts of P as P 3+ or P 5+ , even if the wind is clumped. 
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Wavelength [A 



Fig. 1. — Examples of fits to the Pv wind profiles of 04-09.5 supergiants. Since the red component of the 
doublet is increasingly contaminated by the blend with SilV A1128, the fit to the blue component is given 
greatest weight for spectral types later than ~07. 
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Wavelength [A 



Fig. 2. — Examples of fits to the Pv wind profiles of luminosity class I-V stars with temperature class 
of 07. In all cases, the red component of the doublet is blended with Si IV A1128. The fits to HD 47839 
only provide upper limits to the wind contribution to the Pv wind profiles, which are dominated by the 
underlying photosplreric component. 
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Fig. 3. — Comparison of M with Mq( P 4+ ). The shapes of symbols distinguish radio (circles) and Ha 
measurements, while symbol size separates the primary (large) and secondary (small) samples. Upper limits 
on non-detections are indicated by arrows. Color coding divides the entire sample into early- (02, 03, 03.5), 
mid- (04-07.5), and late-0 types (08-09.7). The dotted line denotes a 1-1 correlation between the two 
measurements . 
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Fig. 4. — q es t = M q(P 4+ )/M(p 2 ) as a function of T e ff for different luminosity classes. Non-detections of P V 
are indicated by downward-pointing arrows, while upward-pointing arrows represent P v detections divided 
by upper limits on non-detected radio or Ha emission. The upper limit at T e g=33 kK in panel for dwarfs is 
for HD 13268, which is a very rapid rotator with a nitrogen-riclr atmosphere. 
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Table 1 

Galactic O-Star Sample 



HD/HDE 


Spectral Type 




Derived Parameters 


v sin i a 


XJ oo 




Notes b 




Classification 


Ref. 


T e ft [kK] 


R/Rq 


Ref. 


km s -1 


km s — 1 


Ref. 




Primary Sample 

13268 ON8 V 


1 


33.0 


10.3 


2 


309 


2150 


3 




14947 


05 If+ 


4 


37.5 


16.8 


2 


133 


2350 


3 




15558 


05 111(f) 


4 


41.0 


18.2 


2 


123 


2800 


3 


SBIO 


24912 


07.5 HI ( n) ((f)) 


4 


35.0 


14.0 


2 


213 


2450 


3 


£ Per 


30614 


09.5 la 


4 


29.0 


32.5 


2 


129 


1550 


3 


a Cam 


34656 


07 11(f) 


5 


34.7 


24.1 


6 


91 


2150 


7 




36861A 


08 III((f)) 


4 


33.6 


15.1 


6 


74 


2400 


3 


A Ori A 


37043A 


09 III 


4 


31.4 


21.6 


6 


116 


2450 


8 


l Ori; SB20 


42088 


06.5 V 


5 


38.6 


10.7 


6 


65 


2200 


3 




47839 


07 V((f)) 


4 


37.5 


9.9 


6 


67 


2300 


8 


SBIO 


66811 


04 I(n)f 


5 


39.0 


19.4 


2 


219 


2250 


9 


C Pup 


93129A 


02 If* 


10 


42.5 


22.5 


2 


180 


3200 


3 


SB 


93250 


03.5 V((f+)) 


10 


46.0 


15.9 


2 


107: 


3250 


3 




149757 


09.5 Vnn 


11 


32.0 


8.9 


2 


372 


1500 


12 


C Oph 


188209 


09.5 lab 


5 


31.0 


19.6 


6 


92 


1650 


7 




190429A 


04 If+ 


4 


39.2 


20.8 


6 


105 


2300 


8 




190864 


06.5 111(f) 


4 


37.0 


12.3 


2 


88 


2500 


3 


SB1? 


191423 


09 III:n 


4 


32.5 


12.9 


2 


436 


1150 


9 




192639 


07 Ib(f) 


5 


35.0 


18.7 


2 


96 


2150 


3 




193514 


07 Ib(f) 


5 


34.5 


19.3 


2 


94 


2200 


3 




193682 


04: 111(f) 


13 


40.0 


13.1 


2 




2800 


9 




203064 


07.5 III:n((f)) 


5 


34.5 


15.7 


2 


305 


2550 


9 


68 Cyg 


207198 


09 Ib-II 


5 


33.0 


16.6 


2 


91 


2150 


3 




209975 


09.5 lb 


5 


32.0 


22.9 


2 


95 


2050 


3 


19 Cep 


210809 


09 lab 


4 


31.5 


21.2 


2 


117 


2100 


3 




210839 


06 I(n)fp 


4 


36.0 


21.1 


2 


219 


2250 


3 


A Cep 


217086 


07 Vn 


4 


36.0 


8.6 


2 


332 


2550 


3 




303308 


04 V((f+)) 


10 


41.0 


11.5 


2 


111 


3100 


3 


SBl? 


Secondary S, 
36486A 


ample 
09.5 II 


5 


30.6 


17.7 


14 


144 


2000 


15 


5 Ori A; SBIOE 


37742 


09.7 lb 


5 


30.5 


22.1 


14 


124 


2100 


8 


C Ori 


46150 


05 V((f)) 


10 


40.9 


11.2 


14 


111 


2900 


8 


SB2? 


46223 


04 V((f+)) 


10 


42.9 


12.4 


14 


82 


2800 


8 




57061 


09 II 


5 


30.6 


17.7 


14 


120 


1960 


15 


r CMa; SBE 


149038 


09.7 lab 


5 


30.5 


22.1 


14 


86 


1750 


8 


H Nor 


149404 


09 la 


4 


31.4 


21.8 


14 


100 


2450 


15 


SB20 


151804 


08 Iaf 


4 


33.2 


21.1 


14 


104 


1600 


8 




152408 


08: Iafpe 


5 


33.2 


21.1 


14 


85 


960 


15 




152424 


OC9.7 la 


5 


30.5 


22.1 


14 


86 


1760 


15 


SBl? 


164794 


04 V((f)) 


4 


42.9 


12.4 


14 


70 


2950 


8 


9 Sgr; SB2? 


188001 


07.5 Iaf 


5 


35.0 


20.5 


14 


93 


1800 


8 


9 Sge; SBl? 



a From Howarth et al. (1997). 

b Binary status: SB1 — single-line spectroscopic binary; SB2 — double-line spectroscopic binary. Suffixes: O - orbit derived; E — eclipsing 
system; ? — possible binary. 

References. — (1) Mathys 1989; (2) Repolust et al. 2004; (3) Haser 1995; (4) Walborn 1973; (5) Walborn 1972; (6) Markova et al. 
2004; (7) Howarth et al. 1997; (8) Groenewegen et al. 1989; (9) Puls et al. 1996; (10) Walborn et al. 2002; (11) Lesh 1968; (12) this work; 
(13) Garmany & Vacca 1991; (14) Martins et al. 2005; (15) Prinja et al. 1990. 
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Table 2 

Log of Far-UV Observations 



HD/HDE 


Spectral Type 


Obs. a 


Data Set b 


PI 


Date Obs. c 


Int. Time d 


N(exp.) 


e Mode f 


Aperture 


Notes 


Primary Sample 

13268 ON8 V 


F 


P1020304 


Jenkins 


1999-11-24 


4438 


2 


TTAG 


LWRS 




14947 


05 If+ 


F 


E0820201 


Massa 


2004-09-30 


5614 


3 


TTAG 


LWRS 




15558 


05 111(f) 


F 


P1170101 


Hutchings 


1999-11-28 


2472 


4 


HIST 


LWRS 




24912 


07.5 III(n)((f)) 


c 


C017-001 




1972-10-30 


36.32 


18 


U2 






30614 


09.5 la 


c 


C021-001 




1972-11-09 


51.59 


29 


U2 






34656 


07 11(f) 


F 


P1011301 


Savage 


2000-03-01 


4179 


7 


HIST 


LWRS 




36861A 


08 III((f)) 


c 


C026-003 




1972-11-27 


79.96 


32 


U2 






37043A 


09 III 


c 


C033-005 




1972-12-25 


70.01 


28 


U2 






42088 


06.5 V 


F 


P1021101 


Jenkins 


2000-11-05 


4219 


9 


HIST 


LWRS 




47839 


07 V((f)) 


c 


C028-004 




1972-12-04 


28.58 


19 


U2 






66811 


04 I(n)f 


c 


C044-001 




1973-02-22 


71.91 


37 


U2 






93129A 


02 If* 


F 


P1170202 


Hutchings 


2000-01-27 


7371 


2 


TTAG 


LWRS 




93250 


03.5 V((f+)) 


F 


P1023801 


Jenkins 


2000-02-04 


4140 


4 


HIST 


LWRS 




149757 


09.5 Vnn 


c 


C002-004 




1972-09-05 


48.51 


42 


U2 






188209 


09.5 lab 


F 


S5231106 


Sahnow 


2004-10-08 


11542 


24 


HIST 


HIRS 


g 


190429A 


04 If+ 


F 


P1028401 


Jenkins 


2000-07-18 


5390 


10 


HIST 


LWRS 




190864 


06.5 111(f) 


F 


E0820501 


Massa 


2004-05-23 


3468 


3 


TTAG 


LWRS 




191423 


09 III:n 


F 


E0821301 


Massa 


2004-06-09 


4795 


4 


TTAG 


LWRS 




192639 


07 Ib(f) 


F 


C1710101 


Nichols 


2002-09-04 


14720 


7 


TTAG 


LWRS 




193514 


07 Ib(f) 


F 


E0820701 


Massa 


2004-06-08 


6478 


2 


TTAG 


LWRS 




193682 


04: 111(f) 


F 


E0820301 


Massa 


2004-06-09 


5818 


5 


TTAG 


LWRS 




203064 


07.5 III:n((f)) 


c 


C152-004 




1979-07-08 


23.29 


13 


U2 






207198 


09 Ib-II 


F 


P1162801 


Snow 


2000-07-23 


13180 


3 


TTAG 


LWRS 




209975 


09.5 lb 


F 


D0140302 


Federman 


2004-07-26 


1294 


3 


HIST 


LWRS 


h 


210809 


09 lab 


F 


P1223103 


Jenkins 


2000-08-08 


10065 


18 


HIST 


LWRS 




210839 


06 I(n)fp 


F 


P1163101 


Snow 


2000-07-22 


6050 


10 


HIST 


LWRS 




217086 


07 Vn 


F 


E0820801 


Massa 


2004-07-25 


4970 


2 


TTAG 


LWRS 




303308 


04 V((f+)) 


F 


P1221602 


Jenkins 


2000-05-27 


7692 


12 


HIST 


LWRS 




Secondary S, 
36486A 


ample 
09.5 II 


c 


C025-001 




1972-11-21 


70.76 


48 


U2 






37742 


09.7 lb 


c 


C024-004 




1972-11-18 


71.94 


44 


U2 






46150 


05 V((f)) 


F 


P1021401 


Jenkins 


2001-03-05 


4888 


9 


HIST 


LWRS 




46223 


04 V((f+)) 


F 


C1680302 


Bruhweiler 


2004-02-23 


7122 


14 


HIST 


LWRS 




57061 


09 II 


c 


C046-001 




1973-03-16 


71.37 


43 


U2 






149038 


09.7 lab 


F 


P3032601 


Williger 


2004-09-11 


58 


1 


HIST 


LWRS 


h 


149404 


09 la 


F 


P1161702 


Snow 


2001-08-07 


17850 


38 


HIST 


LWRS 




151804 


08 Iaf 


c 


C065-002 




1974-09-05 


13.47 


11 


U2 






152408 


08: Iafpe 


c 


C067-002 




1974-09-04 


9.43 


8 


U2 






152424 


OC9.7 la 


F 


Z9016301 


Andersson 


2002-04-18 


7313 


2 


TTAG 


LWRS 




164794 


04 V((f)) 


B 


BEFS1037 




1993-09-16 


257 


1 


TTAG 






188001 


07.5 Iaf 


F 


E0820901 


Massa 


2004-05-29 


58 


1 


HIST 


LWRS 


h 



a Observatory used to obtain spectra: B = BEFS ; C= Copernicus; F—FUSE. 
b The identifier for the data set in MAST. Copernicus data sets are coadded scans. 

C UT date when the scan or integration started, expressed as year-month-day. 

d For Copernicus spectra: the time betwen the start of the first scan and the end of the last scan, in hours. For FUSE and BEFS spectra: 
the total integration time, in s. 

e For Copernicus spectra: number of scans in coadded file. For FUSE amd BEFS spectra: number of exposures comprising the total 
integration. 

f For Copernicus spectra: the detector. For FUSE and BEFS spectra: the mode of detector operation. 

S FP-SPLIT observation; part of a test for observing bright objects. 
h LiFl only observation. 
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Table 3 

Rescaled Measurements of M(radio) 



HD 


Spectral Type 




Previous 








Revised 






dold. 

(kpc) 


log M 
(M 0 /yr) 


Ref. 


M v 


Ref. 


dnew 

(kpc) 


f corr 


log M 
(M 0 /yr) 


Primary Sample 


















14947 


05 If+ 


2.0 


-5.52 


1 


-5.94 


2 


2.0 


1.0 


-5.52 


30614 


09.5 la 


1.0 


-5.41 


3 


-7.00 


2 


1.2 


1.3 


-5.30 


36861A 


08 III((f)) 


0.5 


-6.04 


3 


-5.85 


4 


0.7 


1.6 


< -5.84 


37043A 


09 III 


0.5 


-6.50 


3 


-6.24 


4 


0.6 


1.3 


-6.39 


66811 


04 I(n)f 


0.4 


-5.54 


3 


-6.32 


2 


0.5 


1.4 


-5.39 


149757 


09.5 Vnn 


0.1 


-7.41 


3 


-4.35 


2 


0.2 


2.8 


-6.96 


190429A 


04 If+ 


1.7 


-5.34 


1 


-6.51 


4 


2.2 


1.5 


-5.16 


210839 


06 I(n)fp 


0.8 


-5.65 


3 


-6.40 


2 


0.9 


1.2 


-5.57 


Secondary 


Sample 


















36486A 


09.5 II 


0.5 


-5.97 


3 


-5.73 


5 


0.4 


0.7 


-6.12 


37742 


09.7 lb 


0.5 


-5.60 


3 


-6.28 


5 


0.4 


0.7 


-5.75 


57061 


09 II 


1.5 


-5.20 


3 


-5.77 


5 


0.9 


0.5 


-5.50 


149038 


09.7 lab 


1.3 


-5.36 


3 


-6.28 


5 


1.0 


0.7 


< -5.51 


149404 


09 la 


1.4 


-4.91 


3 


-6.29 


5 


0.9 


0.5 


-5.21 


151804 


08 Iaf 


1.9 


-5.00 


3 


-6.30 


5 


1.2 


0.5 


-5.30 


152408 


08: Iafpe 


1.9 


-4.87 


3 


-6.30 


5 


1.4 


0.6 


-5.09 


152424 


07.5 Iaf 


1.9 


-5.26 


3 


-6.28 


5 


1.4 


0.6 


-5.48 



References. — (1) Scuderi et al. 1998; (2) Repolust et al. 2004; (3) Lamers Sz Leitherer 1993; (4) Markova 
et al. 2004; (5) Martins et al. 2005. 



Table 4 

Properties of the Pv Resonance Doublet 



Property 




Value 


Ref. Remarks 


Rest Wavelength 


^0,blue , ^0,red 


1117.977 A, 1128.008 A 


1 


Oscillator Strength 


fblue, f red. 


0.473, 0.234 


1 


Abundance by number with respect to H 


12.00 + log(N P /N H ) 


5.45 ± 0.06 


2 Solar Photosphere 



References. — (1) Morton 2003; (2) Biemont et al. 1994. 
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Table 5 

Derived and Adopted Mass-Loss Rates 



HD/HDE 


Spectral Type 


0 


logM(j(P 44 


) a 


logM(radio) a 






LogM(Ho;) a 




log Mq 


a 


log M 


a 


Ref. 


logM 


<7 


Ref. 


Primary Sample 






















13268 


ON8 V 


0.7 


-8.14 ( 


-0.30 


+0.18) 








< -6.59 


(-0.34, +0.16) 


1 


14947 


05 If+ 


1.0 


-6.58 ( 


-0.12 


+0.10) 


-5.52 ( 


-0.22, +0.19) b 


2 


-5.07 


(-0.10, +0.10) 


1 


15558 


05 111(f) 


1.0 


-7.64 ( 


-0.30 


+0.18) 








-5.25 


(-0.19, +0.11) 


1 


24912 


07.5 III(n)((f)) 0.8 


-7.21 ( 


-0.12 


+0.10) 








-5.97 


(-0.19, +0.12) 


1 


30614 


09.5 la 


1.0 


-7.43 ( 


-0.12 


+0.10) 


-5.30 


±0.18 


3 


-5.22 


(-0.10, +0.10) 


1 


34656 


07 11(f) 


1.2 


-7.36 ( 


-0.12 


+0.10) 








-5.62 


(-0.09, +0.19) 


4 


36861A 


08 III((f)) 


0.7 


-7.88 ( 


-0.30 


+0.18) 


< -5.84 


±0.18 


3 


-6.01 


(-0.10, +0.24) 


4 


37043A 


09 III 


0.7 


-8.20 ( 


-0.30 


+0.18) 


-6.39 


±0.18 


3 


-5.86 


(-0.10, +0.10) 


4 


42088 


06.5 V 


1.0 < 


-9.34 












-6.21 


(-0.10, +0.08) 


4 


47839 


07 V((f)) 


0.7 < 


-8.79 












-5.92 


(-0.10, +0.10) 


4 


66811 


04 I(n)f 


0.5 


-6.36 ( 


-0.12 


+0.10) 


-5.39 


±0.15 


3 


-5.06 


(-0.13, +0.13) 


1 


93129A 


02 If* 


1.0 


-7.14 ( 


-0.12 


+0.10) 








-4.58 


(-0.10, +0.10) 


1 


93250 


03.5 V((f+)) 


1.0 


-7.42 ( 


-0.12 


+0.10) 








-5.46 


(-0.16, +0.11) 


1 


149757 


09.5 Vnn 


1.0 < 


-9.88 






-6.96 


±0.16 


3 


< -6.75 


(-0.66. +0.17) 


1 


188209 


09.5 lab 


1.0 


-8.04 ( 


-0.30 


+0.18) 








-5.80 


(-0.10, +0.08) 


4 


190429A 


04 If+ 


0.7 


-6.36 ( 


-0.12 


+0.10) 


-5.16 ( 


-0.16, +0.15) b 


2 


-4.85 


(-0.10, +0.08) 


4 


190864 


06.5 111(f) 


0.7 


-7.15 ( 


-0.12 


+0.10) 








-5.86 


(-0.18, +0.16) 


1 


191423 


09 III:n 


0.7 


-7.96 ( 


-0.30 


+0.18) 








< -6.39 


(-0.19, +0.11) 


1 


192639 


07 Ib(f) 


0.7 


-6.75 ( 


-0.12 


+0.10) 








-5.20 


(-0.10, +0.10) 


1 


193514 


07 Ib(f) 


0.7 


-7.00 ( 


-0.12 


+0.10) 








-5.46 


(-0.13, +0.12) 


1 


193682 


04: 111(f) 


0.7 


-7.06 ( 


-0.12 


+0.10) 








-5.76 


(-0.28, +0.14) 


1 


203064 


07.5 III:n((f)) 


1.0 < 


-8.42 












-5.85 


(-0.20, +0.13) 


1 


207198 


09 Ib-II 


0.7 


-7.86 ( 


-0.30 


+0.18) 








-5.75 


(-0.17, +0.15) 


1 


209975 


09.5 lb 


1.0 


-8.37 ( 


-0.30 


+0.18) 








-5.67 


(-0.17, +0.10) 


1 


210809 


09 lab 


0.7 


-7.44 ( 


-0.12 


+0.10) 








-5.28 


(-0.12, +0.12) 


1 


210839 


06 I(n)fp 


0.7 


-6.63 ( 


-0.12 


+0.10) 


-5.57 


±0.15 


3 


-5.16 


(-0.16, +0.16) 


1 


217086 


07 Vn 


1.0 < 


-8.67 












< -6.64 


(-0.68, +0.19) 


1 


303308 


04 V((f+)) 


0.7 


-7.14 ( 


-0.12 


+0.10) 








-5.79 


(-0.13, +0.13) 


1 


Secondary Sample 






















36486A 


09.5 II 


1.0 


-8.32 ( 


-0.30 


+0.18) 


-6.12 


±0.15 


3 








37742 


09.7 lb 


1.0 


-8.26 ( 


-0.30 


+0.18) 


-5.75 


±0.15 


3 








46150 


05 V((f)) 


1.0 < 


-9.21 












< -8.00 




5 


46223 


04 V((f+)) 


1.0 < 


-9.06 












-5.80 




5 


57061 


09 II 


1.0 


-8.05 ( 


-0.30 


+0.18) 


-5.50 


±0.15 


3 








149038 


09.7 lab 


1.0 


-8.05 ( 


-0.30 


+0.18) 


< -5.51 


±0.18 


3 








149404 


09 la 


0.7 


-7.27 ( 


-0.12 


+0.10) 


-5.21 


±0.15 


3 








151804 


08 Iaf 


1.5 


-7.09 ( 


-0.12 


+0.10) 


-5.30 


±0.17 


3 








152408 


08: Iafpe 


1.5 


-6.99 ( 


-0.12 


+0.10) 


-5.09 


±0.15 


3 








152424 


OC9.7 la 


0.7 


-7.59 ( 


-0.30 


+0.18) 


-5.48 


±0.17 


3 








164794 


04 V((f)) 


1.0 < 


-9.15 












-5.68 




5 


188001 


07.5 Iaf 


0.7 


-6.98 ( 


-0.12 


+0.10) 








-5.29 




5 



a In units of M@/yr. 

b The uncertainties quoted by Scuderi et al. 1998 have been increased by ±0.10 dex to include uncertainties in the 
distance. 

References. — (1) Repolust et al. 2004; (2) Scuderi et al. 1998 (revised in Table 3); (3) Lamers Sz Leitherer 1993 (revised 
in Table 3); (4) Markova et al. 2004; (5) Puls et al. 1996 (their Table 11). 
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a Median value of M(Pv) s /M(p 2 ) s for the range of spectral types. 
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